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Abstract 

We address the propagation of twin-beam of radiation through Gaus- 
sian phase-sensitive channels, i.e. noisy channels with squeezed fluctua- 
tions. We find that squeezing the environment always reduces the survival 
time of entanglement in comparison to the case of simple dissipation and 
thermal noise. We also show that the survival time is further reduced if 
the squeezing phase of the fluctuations is different from the twin-beam 
phase. 

1 Introduction 

The use of phase-sensitive environments has been addressed by many authors 
012 El El for preservation of macroscopic quantum coherence. In fact, by adding 
squeezed fluctuations to dissipation, superpositions of states preserve their co- 
herence longer than in presence of dissipation alone. A question arises whether 
or not phase-sensitive environments may be also used to preserve entanglement. 
In this paper, in particular, we study the behavior of a twin-beam state of radi- 
ation (TWB) propagating through a Gaussian noisy phase-sensitive channel, in 
order to investigate whether or not squeezing the environment is useful to pre- 
serve continuous variable entanglement. The answer turns out to be negative. 
The survival time of entanglement in a squeezed bath is always smaller than in 
purely dissipative or thermal ones, and the degradation is more pronounced as 
far as the bath fluctuations are out of the TWB phase. 



2 TWB in a Gaussian bath 

The propagation of a TWB interacting with a general Gaussian environment, 
can be modeled as the coupling of each part of the state with a non-zero tem- 
perature squeezed reservoir. The dynamics can be described by the two-mode 
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Master equation 



^ = {T(l + N)L[a]+T(l + N)L[b]+TNL[at]+TNL[tf] 

+TMM[a)] + TM*M[a] + TMM[tf] + TM*M[b}}p t , (1) 

where p t = p(t) is the system's density matrix at the time t, T is the damping 
rate and N and M are the effective photons number and the squeezing parameter 
of the bath respectively (which are assumed to be equal for the two channels). 
L[0] is the Lindblad superoperator, L[0]p t = Op t O^ — ^O^Opt — and 
M[0]p t — Op t O — ^OOpt — \ptOO. Of course, the dynamics of the two modes 
are independent on each other. 

Using the differential representation of the superoperators in equation 
the corresponding Fokker-Planck equation for the two-mode Wigner function 
W ee W(x 1 ,y 1 ;x 2 ,V2) is given by 



[ 3=1 i,3=l J 



(2) 



where, for the sake of simplicity, we put x = (x±, yi; x 2 , 2/2) = (£1,2:2; 2:3, 2:4), 
t = Tt/j and 7 = (2N + l) -1 . In equation (0) a,j(x) and dij are the matrix 
elements of the drift and diffusion matrices A(x) and D respectively, which are 
given by 

A(i)=~7i?L, (3) 
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Notice that the drift term is linear in x and the diffusion matrix does not depend 
on x. The positivity of D requires that \M\ < (2N + l)/2. Moreover, to ensure 
positivity of the density matrix p t > we need |M| 2 < N(N+1), which includes 
the former condition. 

The solution of the Fokker-Planck © can be calculated analytically. For 
the case 3m[il/] = 0, and considering (without loss of generality) TWB with 
real parameter as the initial state, i.e. po ee ptwb = |TWB))((TWB|, where 
|TWB)) = V 1 ~ £ 2 E„ C p \p)\p), £ € K, the solution assumes the simple form 
□ 
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where £ 2 = T,?(r,T,rhh,n s ), j = 1,2,3,4, are 

Sf = o\ e- rt + D\(t) , £2 = a 2_ e -n + D 2_ (t) ^ 

(6) 

£§ = a 2_ e -rt + D 2 +{t) ; S 2 = a 2 + e -rt + D 2_ {t) i 
with o± = \e ±2X , £ = tanhA, and 

D 2 ± (t) = (l-e~ rt ). (7) 

For the general case (M complex) the analytical solution of J5J is quite cum- 
bersome, and we do not explicitly write it here. 

Notice that if we suppose the environment composed by a set of oscillators 
excited in a squeezed-thermal state of the form v — S(Qpt\ 1 S' ! ((), with £ = 

|C|e 4e , S(C) =exp{i[C* a t2 -Ca 2 ]} and p th = (l + n th )- 1 [n th /(l + n th )] ata , then 
we can rewrite the parameters N and M in terms of the squeezing and thermal 
number of photons n s = sinh 2 |£| and n t h respectively. We have M = \M\e lS [H] 
and 



\M\ = (1 + 2 n t h) y/n s {\ + n s ) and N = n th + n s (l + 2 n th ) • (8) 
This parametrization automatically guarantees the semipositivity of pt ■ 



3 Separability 

A quantum state of a bipartite system is separable if its density operator can be 
written as g = J2kP k(Tk ® Tk i wnere {Pk} is a probability distribution and t's and 
a's are single-system density matrices. If a state is separable the correlations 
between the two systems are of purely classical origin, otherwise it is entangled. 
A necessary and sufficient condition for separability of Gaussian states is the 
positivity of the density matrix g T , obtained by partial transposition of the orig- 
inal density matrix (PPT condition) [3011]. Notice that the Wigner function of 
a twin-beam is Gaussian and the evolution in a Gaussian environment preserves 
such character. Therefore, we are able to characterize the entanglement at any 
time and discuss its degradation as a function of bath's parameters. The PPT 
condition for a density matrix can be rephrased as a condition on the covariance 
matrix V of the two-mode Wigner function W(xi, y\\ x%, y^)- After defining 

and 

V pk = (A£ p A£ fc ) - f A£ p AC fc W(0 , (10) 

ji 4 
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with A£j = £j — and £ = {xi, yi, Xa, 2/2}, we have that a state is separable 
iff 

s = v + -rj>o (11) 

4 - 

For the state 10, the condition l|ll|l rewrites as |S] 

— , e 2 . e? > - 

16 2 3 ~ 16 



E 2 E 2 > - , E 2 E 2 > - . (12) 



We remind that in this case the squeezing parameter M is real i.e. it has the 
same phase of the TWB parameter (6 = 0). By solving inequalities 1)12(1 with 
respect to time t, we find that the TWB becomes separable for t > t s , where 
the survival time t a = i s (A, T, nth, n s ) is given by 

(13) 




where we have defined 

, m v (1 + 2 n th ) [l + 2n th -e- 2A (l + 2n s )] 

/ = /(A,n th ,n s ) = — . (14) 

4n t h(l + nth) 

As one may expect, t s decreases as nth and n s increase. Moreover, in the limit 
n s — > 0, the threshold time reduces to the value to pertaining a non squeezed 
bath p]. In order to see the effect of squeezing the bath on the survival time 
we introduce the function 

G(X,n th ,n s ) = — -• (15) 
to 

G > means that squeezing leads to a longer survival time, shorter otherwise. 
Results are illustrated in Fig. where we plot G as a function of n s for different 
values of n t h and A. Since G is always negative, we conclude that coupling a 
TWB with an in-phase squeezed bath destroys entanglement faster than the 
coupling with a non squeezed environment. 

We have also evaluated the threshold time for separability in case of an 
out-of-phase squeezed bath i.e. for complex M = \M\ e . The analytical ex- 
pression is quite cumbersome and will not be reported here. However, in order 
to investigate the positivity of S as a function of 9, it suffices to consider the 
characteristic polynomial qs (x) associated to S, and study the sign of its roots. 
A numerical analysis shows that this polynomial has four real roots and three 
of them are always positive. We focus our attention on the other one. In Fig. 
[21 we plot the characteristic polynomial for different values of the parameters 
n t h, n s and 9 (we put e _r * = 0.55): it is apparent that by adding thermal noise 
and non classical fluctuations the sign of the smallest root changes from nega- 
tive (entangled state) to positive (separable state). In other words, the survival 
time becomes shorter. Moreover, in Fig. |3|we show that increasing 9 from 
to 7r/2 the threshold time is further reduced. The behavior of the polynomial 
roots for different values of A and Ft is analogue. 



4 



4 Conclusions 



In this paper we have analyzed the propagation of a TWB through Gaussian 
phase-sensitive noisy channels, and have evaluated the threshold (survival) time 
for the state to become separable. 

We found that the survival time in a squeezed environment is always shorter 
than in a purely dissipative or thermal ones. In addition, the survival time is 
further reduced if the squeezing phase of the fluctuations is different from the 
TWB phase. 
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Figure 1: Plots of the ratio G = (t s — to) /to as a function of the number of 
squeezed photons n s for different values of the TWB parameter A and of the 
number of thermal photons nth when M is real. The values of nth are nth = 10~ 3 
(left) and n t h = 1 (right), while the solid lines, from bottom to top, refer to A 
varying between 0.1 to 1.0 with steps of 0.15. 
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Figure 2: Plots of the characteristic polynomial qs(x) associated to S for a 
TWB propagating in a non classical environment with squeezing parameter 
M = \M\ e w . We set e" r * = 0.55, A = 1, and (a) n t h = n s = 0, (b) n th = 0.5, 
n s = 0, (c) n t h = 0.5, n s = 0.07 and 9 — 0, (d) n t h = 0.5, n s = 0.07 and 
9 = tt/5. Notice that there are always four real roots and three of them are 
positive. The inset is a magnification of the region near to 0: the presence of 
thermal noise and non classical fluctuations with non zero phase reduces the 
survival time. 
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Figure 3: Plots of the characteristic polynomial qs{x) associated to S for a 
TWB propagating in a non classical environment with squeezing parameter 
M — \M\ e . This plot only shows the region near to x = (the other three 
roots are always positive). We set e~ rt = 0.55, A = 1, n t h — 0.5, n s = 0.07 and, 
from left to right, 9 — 0, 7r/10, tt/5 and tt/2: a non real squeezing parameter of 
the bath always reduces the survival time. 
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